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ABSTRACT 

The discovery of hypervelocity binary stars (HVBs) in the Galactic halo would provide 
definite evidence of the existence of a massive black hole companion to Sgr A* . Here we 
use an hybrid approach to compute the rate of ejection and the total number of HVBs 
produced by a hypothetical intermediate-mass black hole (IMBH, M2 < 10^ M©) 
orbiting Sgr A* . Depending on the mass of M2 and on the properties of binary stars 
in the central parsec of the Milky Way, we show that the number of undisrupted HVBs 
expected to be expelled from the Galactic Center before binary black hole coalescence 
ranges from zero to a few dozens at most. Therefore, the non-detection of stellar 
binaries in a complete survey of hypervelocity stars would not rule out the occurrence 
of an IMBH-Sgr A* in-spiralling event within the last few x 10* years. 

Key virords: black holes physics - Galaxy: center - stellar dynamics 



1 INTRODUCTION 

Hypervelocity stars (HVSs) are a natural consequence of 
the presence of a massive black hole (MBH) in the Galactic 
Center (GC). At present several HVSs are known to travel 
in the halo of the Milky Way (MW) with Galactic rest- 
frame velocities between -1-400 and 4-750 kms~^ (Brown et 
al. 2005, 2006, 2007). Only the tidal disruption of a tight 
stellar binary by a single MBH in Sgr A* or the scatter- 
ing of a single star by a hypothetical MBH binary (MBHB) 
can kick a 3-4 M© star to such extreme velocities (e.g. Hills 
1988; Yu & Tremaine 2003, hereafter Y03). Direct obser- 
vational evidence for a secondary intermediate-mass hole 
(IMBH) closely orbiting Sgr A* is difficult to establish, how- 
ever. In Sesana et al. (2007) we showed that the observed 
velocity distribution of HVSs appears to marginally disfa- 
vor the MBHB ejection mechanism, though the statistics 
is still rather poor. Lu et al. (2007, hereafter L07) showed 
that tight binary stars can be ejected by a MBHB without 
being tidally torn apart: the discovery of just one hyperve- 
locity binary star (HVB) in forthcoming deep stellar surveys 
could then provide evidence of the existence of a massive or 
intermediate-mass black hole companion to Sgr A*. 

The analysis of L07 is the starting point of this Letter. 
Our goal is to provide an estimate of the number of HVBs 
expected to be produced by the in-spiral of an IMBH onto 
Sgr A* , using the results of scattering experiments between 
a MBHB and a bound stellar cusp discussed in Sesana et al. 
(2008, hereafter SOS). We will show that a short burst of 



HVSs accompanied by a few HVBs would be an incontro- 
vertible signature of a recent in-spiral. By contrast, depend- 
ing on the properties of the population of binary stars in the 
GC, it is possible that a fast binary black hole in-spiral and 
coalescence may occur without the ejection of a single HVB 
in the Galactic halo. 



2 MBHB-STAR INTERACTIONS 

Consider a star of mass m, orbiting the primary hole Mi, 
and assume, for simplicity, that the secondary hole M2 
(m, ^ M2 <C M\) is in a circular orbit of radius a 
around Mi. When the star experiences a close encounter 
with M2, its velocity is of order the MBHB circular veloc- 
ity, V, ^ Vc = {GMi/a)^''^ . A star having closest approach 
distance to M2 equal to rmin.2 ^ a will be subject to a ve- 
locity variation Au» ~ {GM2/rmin,2)^^^ as a result of a (spe- 
cific) force ~ GM2/r'^^^ 2 applied for an encounter timescale 
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(Quinlan 1996; ¥03). This leads to 
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where x = rmin,2/a and q = M2/M1. Since in the limit of 
close energetic encounters the ejection velocity of the star 
is, to first order, Vcj oc An,, equation ([1} shows two im- 
portant scalings: (1), v^j is inversely proportional to the 
square root of the closest approach distance to A/2 dur- 
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ing the interaction; and (2), the closest approach distance 
required to eject a star above a given speed (in units of 
Vc) scales with the MBHB mass ratio q. To verify these 
simple analytical estimates, we have performed 15 sets of 
3-body scattering experiments, using the setup described 
in SOS, for mass ratios q = 1/81, 1/243, 1/729, and eccen- 
tricities e = 0,0.1,0.3,0.6,0.9. In each set we integrated 
5,000 orbits drawn from an isotropic distribution of stars 
bound to Ml, and recorded Ucj, rmin,2, and rmin.i (the clos- 
est approach distance to Mi). The stellar semi-major axis 
a* is randomly sampled from fifty logarithmic bins span- 
ning the range 0.03a < a, < 10a. The stellar specific an- 
gular momentum L, is sampled in the interval [0, L^ nj^xli 
according to an equal probability distribution in Li where 
i-i.tnax = V GM\at is the specific angular momentum of a 
circular orbit of radius a,. A population of stars with such 
distribution in has mean eccentricity (e) — 0.66, corre- 
sponding to an isotropic stellar distribution (e.g. Quinlan, 
Hernquist, & Sigurdsson 1995). We stress that, on average, 
the ejection velocity does not depend on the details of the 
initial Keplerian orbit of the star around Mi, but only on 
rmin,2. Results are plotted in Figure[T]for an assumed MBHB 
eccentricity e — 0.1. The scaling Vc\ oc a;"^''^ breaks down 
for encounters closer than x ~ 0.1 g: the ejection velocity 
tends to (uoj) ^ 3X4, the maximum ejection speed at infin- 
ity predicted by simple arguments on elastic scattering. We 
checked that the precise value of e does not play a significant 
role on determining (uoj). 

The third body in our experiments can be thought of 
either as a single star or as a stellar binary. A binary star 
of mass rrib = m«,i -I- m«,2 and semimajor axis at is broken 
apart by tidal forces if its center-of-mass approaches a com- 
pact object of mass M within the distance (e.g. Miller et al. 
2005) 

/ M\i/^ 

tt ~ 3 ah 

\ vribJ 

- (lo^fc^)'"(aw)- 

Such "breakup" radius must be compared with the closest 
approach distance rmin,2 required for a hypervelocity ejec- 
tion. If TT < »"min,2, a stcUar binary may be kicked to high 
speeds while preserving its integrity (L07). 

Our three-body approximation does not account for the 
internal degrees of freedom of the stellar binary. In partic- 
ular, a strong interaction with the MBHB can result in the 
merger of the two stars. Simulations of stellar binary-binary 
interaction in the context of star cluster dynamics show 
that a merger event is a quite common dynamical outcome 
(Fregeau et al. 2004). However, the dynamical regime we 
consider is different. The stellar binary experiences a com- 
plex weak dynamical interaction with the MBHB (that is 
unlikely to affect the binary internal structure). Ejection or 
breakup is caused instead by an instantaneous strong en- 
counter with one of the two MBHs. Simulations of strong 
three body encounters involving a stellar binary and a sin- 
gle MBH (with parameters similar to those considered here; 
Ginsburg & Loeb 2007) show that the merger of the two 
stars happens at most in ~10% of the cases. 

To make definite predictions, the results of our scat- 
tering experiments must be scaled to the GC. The main 
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Table 1. Parameters of the different models. The quantities 
7, ro, po, 9, O.0, Po, and Vc,o are, respectively, the stellar cusp slope, 
the influence radius of Sgr A* , the stellar density at ro, the MBHB 
mass ratio, its separation when ejections start, the MBHB orbital 
period and circular velocity at ao- 



parameters of our MW models are summarized in Table [T] 
(see Sesana et al. 2007 and SOS for details). The reservoir 
of stars in the central parsec of the MW is well described 
by a power-law density profile, p{r) = po{r /ro)~'^ , around 
a 3.5 X 10® Mq MBH. Here ro is the characteristic radius 
within which the total stellar mass is 2Mi (the "radius of in- 
fiuence" of Sgr A*). As the hypothetical secondary hole M2 
sinks in, it starts ejecting background stars when the total 
stellar mass enclosed in its orbit is Mt(< a) ~ M2 (Mat- 
subayashi et al. 2007). FoUowing SOS, we set the MBHB at 
initial separation ao such that Mt{< ao) = 2Af2. From 7, 
po, ro, and q we can derive the parameters ao, Vc,o, and the 
period P(ao) = Po. We take a velocity threshold for escaping 
the MW potential of Vcsc = 850 kms"^ at ro (e.g. Smith et 
al. 2007). The horizontal lines in Figure [T] depict the quan- 
tity Vbsc/Vc at orbital separation a = ao and a = ao/10, 
while vertical lines mark the tidal disruption radius rT,2 in 
units of a for the same two MBHB separations. An equal- 
mass stellar binary with rrib — 2 Mq, and ai, — 0.1 AU 
was assumed. The figure shows that, on average, stars must 
approach M2 within a distance x < rT,2/a in order to be 
ejected. Most binary stars will then be tidally disrupted dur- 
ing the strong interaction with M2 , and only a few tight bi- 
naries with af,<0.1 AU may be survive intact and become 
HVBs. Note that, while stellar binaries can be also broken 
apart by Mi (rT,i S> rT,2), it is the secondary hole M2 that 
is largely responsible for dissociating candidate HVBs. This 
is because there is no connection between tidal dissociation 
by Ml and hypervelocity kicks, while a close approach to 
M2, required to gain hypervelocity, can break up the binary. 



3 HYPERVELOCITY STELLAR BINARIES 

To quantify the fraction of binary stars that are not dis- 
rupted by M2 (and Mi), we need to specify their mass and 
semi-major distributions. In our default model, we assume a 
log-flat distribution of semi-major axis, 

p{at)dat = dat/at, (3) 

in the range 10-2 < at < 1 AU (Heacox 1998). The lower 
limit is set by the contact separation of two solar-mass stars, 
while the upper limit considers that binaries with at, > lAU 
are unlikely to survive the dense stellar environment of the 
GC (e.g. Y03). We have also run a case with the distribution 
of semi-major axis arising from a log-normal distribution of 
binary periods Pt: 
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Figure 1. Thick curves: mean stellar ejection velocity (in unit 
of Vc) as a function of the dimensionless minimum distance of 
approach to M2, x = rmin 2/1- The MBHB has an eccentricity of 
e = 0.1 and a mass ratio q = 1/81 {solid line), q = 1/243 {long- 
dashed line) and q = 1/729 (short- dashed line). The horizontal 
lines mark the escape velocity Veac = 850 kms~^ in units of Vc{a) 
for a = ao, a = 0.1 ao, and different mass ratios (using the same 
line styles as above). Similarly, the vertical lines mark the tidal 
disruption ratio 2/c' for an equal-mass binary with = 2 Mq 
and ai, = 0.1 AU (leftmost three lines for a = ao, rightmost 
three for a = 0.1 ao). Dots mark the intersection of corresponding 
horizontal and vertical lines dividing the {vej) /Vc — x plane into 
four quadrants. To produce an HVB, a thick curve must lie in 
the corresponding upper right quadrant, where t)ej > fesc and 
'"miii,2 > 1'T,2- Note how, on average, these stellar binaries tend 
to be disrupted and do not become HVBs. 



p(logPb) = Cexp 



-(logn - {iogn))= 



(4) 



Here C = 0.18 is a normalization constant, (logPb) = 
4.8, crfogFj = 2.3, and the period Pi, is measured in days 
(Duquennoy & Major 1991). The above semi-major axis dis- 
tributions are coupled to two different choices of the stellar 
binary member's mass function (for a total of 4 different 
models). We either assume all binaries to be composed by 
two equal solar-mass stars {rrib = 2M0, default model), or 
to follow a Salpeter initial mass function (IMF) in the range 
1 M0 — 15 Mq for m,.i while m,,2 is randomly chosen in the 
mass range 1 M0 — m,,i. We shall discuss later the effect of 
these different assumptions on our results. 

To estimate the fraction of stellar binaries that survive 
the interaction with the binary black hole and are ejected 
intact as HVBs, we proceed as follows. For fixed q and e, 
we consider orbital separations in the range 0.1 — 1 ao, select 
from our 5,000 simulated orbits those resulting in an ejection 
with Vcj > Vcsc, and denote their number with Nei{a). We 
then assume that each of these "ejection orbits" is followed 
by a binary stellar system with parameters (a;,, m^,i, m*,2) 
drawn from the distributions described above, and calcu- 



Figure 2. Upper panel: fraction of binaries that survive the inter- 
action with the MBHB and are ejected intact with v^j > Vcsc, as 
a function of MBHB separation. Open squares: q = 1/81. Filled 
circles: q = 1/243. Stars: q = 1/729. Lower panel: mean semi- 
major axis ai, of undisrupted HVBs. Equal-mass binaries with a 
log-fiat distribution in at {default model) have been assumed, and 
the MBHB eccentricity is set to e = 0.6. 



late the radii rT,2 and rT,i using equation ([2)). Finally, 
if during the chaotic interaction with the MBH pair it is 
ruiin,2 < rT,2 OT rmin.i < »"T,i, the Stellar binary is counted 
as "disrupted before ejection", and added to 7VTD(a). The 
fraction of HVBs as a function of a is then 



/hvb 



Afq(a) -7VTD(a) 



(5) 



Results are shown in Figure [2] for our default model with 
MBHB eccentricity e = 0.6. The fraction of undisrupted 
HVBs is of order 20 - 40%, dropping to 5 - 20% if the dis- 
tribution of of semi-major axis is derived from equation UJ 
similar fractions are obtained for all the eccentricity values 
we sampled. Surviving hypervelocity binaries have (a6)<0.1 
AU, i.e. only tight binary stars can be ejected undisrupted. 
It is clear from the figure that /hvb and (at) do not signifi- 
cantly change as the MBHB shrinks. We can understand this 
result by noting that {vcj)/Vc oc ^a/rmin,2, i.e. Vcj oc T-~f/_2- 
The ejection velocity (in physical units) does not depend 
then on MBHB separation, but only on the minimum ap- 
proach distance to M2. Figure [2] also shows that the quan- 
tities /hvb and {at) decrease slightly with decreasing black 
hole mass ratios q. This occurs because rT,2 oc Mj^'^ oc q^^^, 
while rmin,2 oc q, i.e. the more massive the secondary hole 
the weaker the interaction required to kick a star above a 
given speed. If q is (say) three times smaller, a binary star 
must approach M2 at a distance three times smaller to be 
ejected. But as the breakup radius rT,2 decreases by just a 
factor 3^''^, fewer tighter stellar binaries can survive undis- 
rupted the tidal field of M2 . 
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4 EJECTION RATES AND DETECTABILITY 

We can now estimate the rate at which binary stars would be 
ejected into the MW halo by an IMBH spiralling into Sgr A* . 
In SOS, we self-consistently computed the orbital evolution 
of such an IMBH in terms of a{t) and e{t) (see figure 8 in 
SOS), and estimated the stellar mass ejection rate dmej/dt. 
Here, we assume that a fraction of scattered stars are 
binaries, and account for the evolving MBHB eccentricity 
during orbital decay by linearly interpolating the fraction 
/HVB(t[, e) along the correct e(a) curve. The ejection rate of 
HVBs can be written as 

1 dmej fb r ff.-. 

(m«) at 2 

where (m,) is the mean stellar mass. Results are shown in 
Figure O for our defauh model and /t = 0.1. The HVB 
ejection rate peaks between 5 x 10~^ — 2 x 10~^ yr~^ over a 
timescale of 10^ — 10^ yr, depending on q. For comparison, we 
also plot the ejection rate of HVSs by the in-spiralling IMBH 
(SOS), as well as the rate of HVSs produced by the tidal 
disruption of a tight stellar binary by a single MBH in Sgr 
A* (Hills' mechanism), as estimated by Y03. In all the cases 
studied, the total number of HVBs, A^Hve is small compared 
to the expected number of HVSs. We find A'hvb = 2S, 9, 4 
for q = 1/81, 1/243, 1/729, respectively. If the distribution 
is log-normal (equation [Jjl , the fraction of tight binaries is 
reduced and the number of HVBs drops by about a factor of 
2. Moreover, in the case of a Salpeter IMF, (m,) ~ 2.5 Mq 
and i?HVB is further reduced by the same factor (see eq. 
[6|. The number of HVBs would trivially increase linearly 
with /(,. The number of ejected hypervelocity binaries is well 
approximated by 

;Vhvb« 280 (170)/,<.^^^J^, (7) 

where /b,<i is the fraction of stars in binaries with < \ 
AU, and 280 (170) is the normalization constant appropriate 
for a log-flat (log-normal) distribution. 

It should be noted that our approach does not account 
for the binary stars that are not initially bound to the 
MBHB and populate its loss cone because of two-body re- 
laxation processes. L07 estimated an HVB ejection rate for 
such unbound population of fewxlO"^ yr~^ in the case of 
a MBHB with q = 0.01 and a = 0.0005 pc. Such a pair 
is expected to have a large eccentricity (e.g. Matsubayashi 
et al. 2007) and a coalescence timescale of only ~ 10^ yr. 
For larger orbital separations, the loss cone is larger but 
the mean ejection velocity is accordingly smaller, leading to 
lower ejection rates. Such rates are one dex smaller than 
those we derived for bound stars and q = 1/81. 



5 SUMMARY 

We have applied the hybrid approach described in SOS to 
compute the rate of ejection and the total number of hy- 
pervelocity binary stars produced by a hypothetical IMBH 
orbiting Sgr A*. Depending on the mass of M2 and on the 
properties of binary stars in the central parsec of the Milky 
Way, we have shown that the number of undisrupted HVBs 
expelled before coalescence ranges from zero to a few dozens 
at most. In particular, we have found that the rapid in-spiral 
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time [yr] 

Figure 3. Ejection rates of HVSs calculated in SOS (thin curves) 
versus the ejection rates of HVBs (thick curves) calculated in 
this work for the default model and = 0.1. Solid lines: Ci = 
0.1. Long-dashed lines: Ci = 0.5. Short-dashed lines: Ci = 0.9. 
The three sets of curves refer, as labeled, to mass ratios q = 
1/81, 1/243, 1/729. The dotted horizontal line marks the ejection 
rate predicted by Hills' mechanism (Y03). 



of a 5 X 10* M© IMBH would generate ~ 40 HVBs, assumme 
a stellar binary fraction of 0.1, m* = 1 Mq, and a log-fiat 
distribution of stellar semi- major axis at- A 10% binary stel- 
lar fraction with <few AU is suggested by numerical sim- 
ulations of dense stellar clusters (e.g Shara & Hurley 2006, 
Portegies-Zwart, McMillan & Makino 2007). The number of 
HVBs is proportional to the mass of the IMBH and inversely 
proportional to m, , so in the case of a top-heavy stellar mass 
function (Schodel et al. 2007), the expected number of HVBs 
would be lower. Moreover, if the distribution of stellar binary 
semi- major axis is log- normal instead of log-fiat, the number 
of tight binaries that can survive a strong interaction with 
the MBH pair is smaller. The combination of these factors 
can potentially decrease the number of expected HVBs to 
zero. 

To conclude, while the observation of even a single HVB 
in the Galactic halo would be a decisive proof of the recent 
in-spiralling of an IMBH into Sgr A*, it is likely that such 
an event would give origin to at most a handful of HVBs. 
Therefore, the non-detection of stellar binaries in a complete 
survey of hypervelocity stars may not be used to rule out 
the existence of an IMBH-Sgr A* pair in the GC. 
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